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Abstract 

We report the first observation of 'orbital truncation rods' — the scattering arising from the 
termination of bulk orbital order at the surface of a crystal. The x-ray measurements, performed 
on a cleaved, single-layered perovskite, Lao.sSri.sMnO^ reveal that while the crystallographic 
surface is atomically smooth, the orbital 'surface' is much rougher, with an r.m.s. deviation from 
the average 'surface' of ~7A. The temperature dependence of this scattering shows evidence of a 
surface-induced second order transition. 
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An issue destined to be of increasing importance in the science of correlated-electrons 
systems in the coming years is the question of whether the electronic behaviour, such as 
charge, orbital or spin order, is different at a surface or interface compared to the bulk of 
the material. This question is of central importance from an applied perspective — any 
device one fabricates from these materials will require electron transport across an interface; 
from a nanoscience perspective — nanoscaled objects are defined by their very large surface- 
to- volume ratio; and from a fundamental perspective — understanding the role of reduced 
dimensionality in determining the electronic behaviour of such systems is a key component 
in understanding the general problem of strongly correlated electron systems. 

In many ways manganites represent an ideal system to address the role of the surface, 
because all the relevant degrees of freedom — charge, spin and orbital — play an active role 
in determining the ground state. As a result, they are exquisitely sensitive to perturbations, 
and hence, might be expected to exhibit relatively large surface effects. In fact, while 
the vast majority of experimental work has focused on bulk properties [l|, Q], the very few 
surface sensitive studies do suggest that the surface of a manganite may be quite different 
from the bulk. For example, the measurements of Freeland et al. showed that the top 
bilayer of a cleaved surface of a ferromagnetic bilayered manganite is magnetically dead, 
though the second bilayer has full bulk magnetic order. Similarly, the half-doped single- 
layer manganite, which is known to exhibit bulk orbital order below Tqo — 231 K, showed 
no superlattice reflection corresponding to such ordering down to 80 K, when probed with 
surface-sensitive LEED |4]. Thus, there is evidence that the surface of manganite crystals 
may exhibit different electronic phases than the bulk — a phenomena that has been coined 
'electronic reconstruction' in analogy with more traditional surface science studies of 
atomic reconstruction driven by surface energetics. 

In this paper, we present the first study of orbital order at a surface, utilizing surface x- 
ray scattering. The advantage of using x-ray scattering in this endeavour is that it provides 
information on both the electronic surface (specifically, atomic displacements associated 
with orbital order) and the chemical surface, separately. It is well known that the abrupt 
truncation of a sample by its surface gives rise to rod-shaped x-ray scattering around the 
Bragg reflections, the so-called crystal truncation rods (CTR). The intensity profiles of such 
rods provide detailed, quantitative information about the chemical surface [7]. Similarly then, 
the abrupt termination of orbital order by a surface should give rise to orbital truncation rods 
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(OTR) running through the orbital superlattice peaks. These would allow one to quantify 
the surface orbital order. Here, we present the first observation of such orbital truncation 
rod scattering. This allows us to quantitatively probe the orbital order in the vicinity of the 
surface. Working with a freshly cleaved LacsSr^MnC^, we show that the orbital order has 
a significantly rougher interface than the chemical surface and a shorter correlation length 
than in the bulk. These results imply that the surface order parameter is reduced relative 
to the bulk. Further, the temperature dependence of the surface orbital order suggests that 
the surface exhibits a more continuous transition than the bulk ordering transition. Taken 
together, these results demonstrate that the surface is a significant perturbation in regard 
to the electronic phenomena of orbital ordering and point to the need for an improved 
theoretical understanding of how the surface couples to such electronic orders. 

We first discuss surface preparation and characterization. Single crystal samples of 
La .5Sr 15 MnO4, which has the layered perovskite structure shown in Fig. [1] (a), were grown 
by the floating zone method. The samples were then either cleaved perpendicular to the 
c-axis in He and placed in a sealed He-filled can, or cleaved in air, and put into a closed-cycle 
refrigerator. Both treatments gave the same results, indicating that the surface is stable to 
brief exposure to air, a result also found elsewhere Figure [D(b) shows an atomic force 
microscope (AFM) picture of a freshly cleaved sample. The terraces seen here are on the 
order of 1 /iin in size, much larger than the coherence length of the x-rays. The step height 
observed here was 6.2A± 0.2 A, corresponding to c/2. We note that similar AFM pictures 
were obtained at the end of a week long x-ray experiment in vacuum. Analysis of the CTR 
intensity allows a determination of the precise termination of the cleaved sample 9| . Figure 
[1] (c) shows the CTR scattering intensity distribution along the (00L) and (HL) directions 
together with the calculated intensity distribution for a La-terminated surface and a Mn- 
terminated surface (FigJT] (a)). The experimental data are well reproduced by assuming a 
predominantly La-terminated surface with very little surface relaxation. For comparison, 
there are conflicting reports in the literature for the termination of thin films, with both La 
and Mn termination reported 10l |. 



As discussed above, the structure of the chemical surface is reflected in the scattering 
along the (hkL) direction, where h, k are integer multiples of reciprocal lattice vectors, that 
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FIG. 1: (a) The crystal structure of Lao.sSri.sMnO^. Green, red and white spheres denote La/Sr, 
Mn, and O atoms, respectively, (b) AFM image of a cleaved sample used in this study. The step 
height was 6.2 ± 0.2 A. (c) (00L) and (11L) CTR scattering intensity distribution together with 
the calculated intensity for La and Mn terminated surfaces. The result of the curve fitting shows 
the chemical surface is predominantly La terminated. 

is rods along the surface normal c*-direction through the allowed, bulk Bragg reflections. 
Similarly, scattering originating from the surface of the orbital ordering must be present 
in 'orbital truncation rods' — rods of scattering along the same direction but through the 
superlattice reflections related to the orbital order, i.e. (h + |, k + j, L) where Q, |, 0) is 



the orbital ordering wave- vector ll[ . One of the primary goals of this work was to observe 



such scattering. To this end, we measured the (20L)±(^0) intensity distribution below 
Tqo =231 K, the bulk orbital order transition temperature. We discuss first transverse 
scans taken across such 'orbital rods' (figure L2ta)). For reference a transverse (2 5k 0.2) scan 
profile across the truncation rod arising from the crystallographic surface is also shown (blue 
data). It is as sharp as the instrumental resolution. As discussed below, this is indicative of 
a smooth atomic surface. The open circles show the profile of the (|t2) bulk orbital order 
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Bragg peak. The finite width observed here shows that the correlation length of the bulk 
orbital order is finite, and in fact, we find £oo k = 410A± 8A, in agreement with earlier 



work 



111 ]. The transverse scan taken at (|t0.2), that is across the expected OTR is shown 



in red. A peak is clearly observed at this position. However, this peak is not, by itself, 
evidence of surface scattering. Such a peak could arise from stacking faults or other defects, 
in the bulk orbital order which would result in a finite correlation length for the orbital order 
along the c-direction and hence a tail along the c*-direction for superlattice reflections. 

We can, however, distinguish surface scattering from such bulk scattering by measuring 
the scattering intensity as a function of incidence angle a. Surface scattered intensity will 
be enhanced at a = 8c, the critical angle for total reflection, as a result of the interference 
between the incident and the scattered waves of the photon 12j . The inset to Fig. Et^a) shows 
the scattered intensity as a function of a at the bulk orbital superlattice reflection (|^2) and 
on the nominal OTR at (||0.3). The intensity at (||0.3) has a maximum around 6c = 0.4° 
while a smooth increase is seen at (\\2) with increasing a, a result of the fact that more of 
the beam is intercepted by the sample at higher angles. Since the exit angle (5 for (||0.3) 
is 3 to 1.5 degrees in this scan region, the observed intensity maximum cannot be due to 
a maximum in the sample illumination and detector acceptance, which happens at a. — [3 
in our experimental configuration. Therefore, we conclude that the scattering observed at 
(2 0.2)±(i|0), shown in Fig. [2] (a) does indeed come from surface scattering and not 
from stacking faults or other bulk effects. This is the first such observation of the orbital 
truncation rod scattering. 

Having established that such scattering originates from the truncation of the orbital 
order, we next discuss what can be learned from this new scattering. As shown in Fig. [2] 
(a), the transverse scan profile of the rod intensity for the chemical surface and for the 
orbital surface are quite different. This difference in the scattering profiles directly reflects a 
difference in the surface roughness. [7] The sharp single peak observed for the chemical CTR 
scattering implies a very flat surface, and is consistent with the atomically smooth chemical 
surface observed in the AFM picture (Fig. [T] (b)). In contrast, the transverse profile of the 
orbital truncation rod exhibits two components; a 'sharp' one and a 'broad' one. Such a line 
shape is expected when a surface — in this case the orbital surface — has a well-defined 

n 

average position (a flat plane) but local, finite deviations from this average[7j]. Thus, without 
recourse to any model-dependent fitting we can immediately conclude from the data of FigfS] 
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FIG. 2: (a) Transverse (k) scans through (2 0.2), (|j2), and (||0.2) reciprocal lattice positions 
shown with blue circles, open circles, and red circles, respectively. The red curve shows the result of 
a two-Gaussian fitting for (||0.2) profile, (inset) The incidence angle dependence of the scattered 
intensity at (tt2) and (tt0.3). The enhancement of the latter at a around 6c = 0.4° shows that 
it results from surface scattering, (b) (20L) and {\\L) truncation rod intensity for samples A and 
B at 100K, 170K and 200K. The solid line shows the calculated intensity distribution given by the 
measured I(20L) for sample A multiplied by exp(— cr 2 q^) with surface roughness parameter a = 6.7 
A (see text). 

(a) that there exists a flat chemical surface and a rougher orbital order surface. 

A rougher surface also has the consequence that the intensity along the truncation rod 
must fall more sharply as a function of L than for a less rough surface. Figure [2] (b) shows the 
intensity of the (20L) and (jjL) rods as a function of L measured on two samples (cleaved 
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(b) 




in He, sample A, and in air, sample B) at 100K, 170K, and 200K. For the orbital rod, the 
intensity of the sharp component integrated in the transverse direction is shown. Note, the 
L-dependence of the (20L) intensity of sample A is slightly different from that of sample B, 
suggesting a slightly different degree of flatness of the chemical surface for the two samples. 
On the other hand, the (§|L) intensity distribution for both samples is very similar and in 
both cases, the slope is steeper than that along (20L). This is consistent with the result in 
Fig. [2] (a): the orbital order surface is rougher than the chemical surface. 

We next quantify this result by use of a model for a nearly flat surface due to Andrews and 
Cowley [3]. In this model, the surface is characterized by two parameters: a, the standard 
deviation of the surface height, and £ sur , the correlation length of the surface (Fig. [3]). The 
truncation rod intensity distribution I(Q) from such a surface is given by: 

I(Q) = I s (Q)+I b (Q), (1) 
with I S (Q) = ^Mlex V (-a 2 q 2 z )S(h,k), (2) 



a*q 

'-(r ur ) 2 [i - 

x exp 



(r r ) 2 [l-exp(-a^)] ( , 2 | fc2 



Aa2q 2 v ■ ■■ , , (3) 

where I s and L, are the intensities for the sharp component and the broad component 
respectively, F(Q) is the structure factor at scattering vector Q, and q z is the c*-component 
of the reduced wavevector (e.g. this is 2nL/c for |L| < 0.5). In order to compare our data to 
the predictions of this model, we replace the delta-function in equation (jSJ) with a Gaussian, 
to account for the finite bulk orbital correlation length, and fit the transverse scans to a 
two-Gaussian lineshape as shown by the red solid curve in Fig. [2](a). The L-dependence of 
the parameters derived from these fits are shown in Figure E(a), which plots the ratio of 
the intensity of the broad component and the sharp component, and in Figure [3](b), which 
shows the widths of the two components, for the 170 K data set. The intensity ratio h/I s 
increases rapidly with L. This ratio depends mainly on a, and hence a may be obtained 
from the plot. The solid line in Fig. E](a) shows the expected behaviour for a = 6.7A, which 
reproduces the data quite well. Note, this same value of o is also consistent with the results 
shown in Fig. [2] (b); The solid line in this latter figure is the calculated intensity distribution 
given by taking the product of the measured J(20L) and exp(— u 2 ql), with <r=6.7A; here, we 
assumed the chemical surface is perfectly flat. The in-plane correlation length £ sur is given 
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FIG. 3: (a) The intensity ratio of the broad component observed in the transverse orbital truncation 
rod scans to the sharp component, as a function of L. (inset) Schematic view of the orbital order 
region around the surface, (b) Peak widths of broad- and sharp- component. Calculated values for 
the model discussed in the text are shown by solid lines in both figures. 

by the transverse peak width and is shown in Fig. [3] (b) as a function of L, together with the 
Independence predicted by eq. (JBJ) for a particular value of £ sur . Using the value of a = 6.7A, 
£ sur is found to be 120A. That is, the surface correlation length is much shorter than that of 
the bulk orbital order = 410A. This suggests that some aspect of the surface disrupts 
the orbital order. This may be due to differing surface energetics resulting from the reduced 
coordination, or perhaps varying Sr concentrations at the site of the cleave. 

We next compare the present results to that of previous reports. LEED measurements on 
the same material Q] showed no superlattice intensity down to 80K, implying little orbital 
order is established in the topmost layer. Our result is consistent with this study. The 
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FIG. 4: Temperature dependence of the (||L)-sharp component intensity for L =2, 0.2, 0.3, and 
0.4. 

surface roughness in the orbital order observed here necessarily requires a small order pa- 
rameter or a small ordered volume fraction at the surface. The same authors also performed 
STM measurements on the same system and found unidentified 'electronic roughness at the 
surface' having an in-plane correlation length of several nanometres and surface roughness 
of 6A. Comparing to the present results, one notes that the length scale of the roughness 
observed in the STM study is very similar to that of the orbital order interface observed in 
this study, and it is tempting to surmise that this 'electronic roughness' may be related in 
some way to the rough orbital interface observed here. However, it is difficult to draw any 
firm conclusions on this point without further study. 

Finally, we briefly discuss the temperature dependence of the orbital truncation rod scat- 
tering. Figure H] shows the temperature dependence of the (|jL)-sharp component intensity 
for L =2, 0.2, 0.3, and 0.4 r.l.u. The temperature dependence of the surface-sensitive OTR 
scattering intensity (L=0.2, 0.3 and 0.4 data) is different from that of the bulk-sensitive 
superlattice reflection. The bulk OO transition was reported to be 'probably second order' 
with a relatively discontinuous transition 13j. In contrast, the temperature dependence of 
the OTR scattering intensity is seen to evolve much more continuously through the transi- 
tion. The bulk superlattice reflection intensity is proportional to the volume of the orbitally 
ordered phase, and to the square of the order parameter, or the Jahn- Teller distortion caused 
by the orbital ordering. In contrast, the OTR intensity is proportional to the area of the 
orbital ordered flat surface and to the square of the order parameter. Since the surface 
roughness of the OO phase at 100K is the same as that at 200 K, then the change in in- 
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tensity must be attributed to a change in the magnitude of the atomic displacement, i.e., 
in the orbital order parameter at the surface. Thus, these data suggest that the surface 
orbital order parameter decays more continuously than the bulk. Such phenomena may be 
related to surface-induced second order phase transitions that have been observed in other 
m ate r ia,sQQ. Future wo rk will investigate this possibiuty iu detai,. 

In conclusion, we have performed a surface x-ray scattering study of a cleaved, single- 
layered, half-doped manganite. The chemical surface was found to be extremely flat with a 
La/Sr layer termination. Importantly, we have observed the scattering arising from trunca- 
tion of the orbital order at a surface for the first time. We find that the orbital surface is 
rougher than the chemical surface, and that its temperature dependence is different from the 
bulk orbital order, indicating a surface-induced second order transition. These results are 
further evidence of the importance of surface effects in determining electronic ground states 
in strongly correlated systems and provide a direct route to probe such effects quantitatively. 
Detailed understanding will require theoretical models capable of dealing with such surface 
effects explicitly. 

Methods 

X-ray diffraction measurements were performed at the beamlines X22C and X21 at the 
National Synchrotron Light Source, and at 6ID at the Advanced Photon Source. The mea- 
surements at X22C and 6ID were performed with 6.5 keV x-rays using a six-circle diffrac- 
tometer for surface scattering experiments and those at X21 were performed with 8.8 keV 
x-rays and a standard four-circle diffractometer. Most of the experiments using the six-circle 
diffractometers were conducted with a fixed-a mode in order to keep the illuminated sample 
volume constant, where a denotes the incidence angle, i.e., the angle between the crystal 
surface and the incident x-ray. The angle a was fixed to 0.6°, which corresponds to a pene- 
tration depth of 600 A, unless otherwise noted. This angle is higher than the critical angle 
for total reflection (9c = OA ). The experiments at X21 were performed in a a = f3 mode, 
where f3 denotes the exit angle. Both modes give very similar results. The temperature was 
controlled by closed-cycle He refrigerators. The sample was mounted either in thermally- 
insulating vacuum or in a He exchange gas. The curve fitting for CTR profile was made with 
the ROD program [l^] with Mo for La/Sr site because the form factor of Mo approximates 
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that of the Lao.25Sro.75 mixture well. 
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